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Abstract 

The generation of a kination-dominated phase by a quintessential exponen¬ 
tial model is investigated and the parameters of the model are restricted so 
that a number of observational constraints (originating from nucleosynthesis, 
the present acceleration of the universe and the dark-energy-density parameter) 
are satisfied. The decoupling of a thermal cold dark matter particle during the 
period of kination is analyzed, the relic density is calculated both numerically 
and semi-analytically and the results are compared with each other. It is argued 
that the enhancement, with respect to the standard paradigm, of the cold dark 
matter abundance can be expressed as a function of the quintessential density 
parameter at the onset of nucleosynthesis. We find that values of the latter 
quantity close to its upper bound require the thermal-averaged cross section 
times the velocity of the cold relic to be almost three orders of magnitude larger 
than this needed in the standard scenario so as compatibility with the cold dark 
matter constraint is achieved. 
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1. Introduction 

A plethora of recent data in 12] has indicated that the two major components of the 
present universe are the Cold (mainly I2j) Dark Matter (CDM) and the Dark Energy (DE) 
with density parameters |T], respectively: 

(a) D C dm = 0.24 ±0.1 and (b) D de = 0.73 ± 0.12, (1.1) 

at 95% confidence level (c.l.). The identification of these two unknown substances consists one 
of the most tantalizing enigmas of the modern cosmo-particle theories. 

As regards CDM, the most natural candidates [I| are the weekly interacting massive 
particles, y’s. The most popular of these is the lightest supersymmetric (SUSY) particle (LSP) 
‘5;. However, the extra dimensional (ED) theories give rise to new CDM candidates [6, ,2E]- 
According to the standard cosmological scenario (SC) 0, x’s (i) are produced through thermal 
scatterings in the plasma, (ii) reach chemical equilibrium with plasma and (iii) decouple from 
the cosmic fluid during the radiation-dominated (RD) era (note that these assumptions are, 
also, naturally valid in the case of the so-called second Randall-Sundrum m model, provided 
that the brane-tension is constrained to rather high values El)- The viability of other CDM 
candidates (like axions JT2], axino m , gravitino M, quintessino m) requires a somehow 
different cosmological set-up, which we do not consider in our analysis. In light of eq. CCD, 
the y-relic density, D x /ig, is to satisfy the following range of values: 

(a) 0.09 < n x h 1 2 0 and (b) Q x h 2 0 < 0.13. (1.2) 

As regards DE, quintessence m, a slowly evolving scalar field, has recently attracted 
much attention (for reviews, see ref. |E]). The scalar field is supposed to roll down its potential 
undergoing three phases during its cosmological evolution: Initially its kinetic energy, which 
decreases faster than the radiation, dominates and gives rise to a possible novel period in 
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the universal history termed “kination” CHI. Then, the scalar held freezes to a value close 
to Planck scale and by now its potential energy, adjusted so that eq. (fTTb l is met, becomes 
dominant. Such an adjustment, which certainly does not resolve satisfactorily the coincidence 
problem, is unavoidable in quintessential models (for related suggestions, see refs. HI KUlh 
Other shortcomings such as the lightness of the scalar held m or the time variation of the 
gauge coupling constants |22] are currently under investigation. 

Be that as it may, the viability of a quintessential scenario can be controlled by imposing 
some observational constraints ESI, arising from nucleosynthesis, acceleration of the universe 
and the DE density parameter. Unfortunately no full-satisfactory potential exists, to date 
(for comparative explorations of various potentials, see refs. 1241I25| 1. E.g., the inverse power 
potential m although provides a tracker-type solution m does not ht well m the present- 
day value mm of the quintessence-equation-of-state parameter. Phenomenologically more 
robust ED is the supergravity-inspired EH OH potential without, however, to allow a zero 
minimum of the potential EEIE]. Also, in both cases, the generation in the early universe 
of a kination-dominated (KD) expansion consistent with the fulfillment of the requirements 
above is rather questionable [SOI ED- For these reasons, we decide to examine the simplest 
exponential potential ,22102], which, although does not possess a tracker-type solution mm- 
it can produce a viable present-day cosmology in conjunction with the domination of an early 
KD era, for a reasonable region of initial conditions Esmacsi. 

The departure from the SC, caused by the implementation of a quintessential KD epoch 
can modify the D x /ig calculation, which (as, already, emphasized ei ms m mm) crucially 
dependents on the adopted assumptions. If the quintessential KD phase dominates over 
the radiation (a condition indispensable for the quintessential inflationary model-building 
EH, 2211101), the assumption (iii) of the SC is lifted (note that the assumptions (i) and (ii) are 
maintained). As a consequence, an increase to D x Aq with respect to (w.r.t) its value in the 
SC is implied. This phenomenon was first pointed out in ref. m and was explored in ref. m 
for the parameters of the exponential potential, which support a global attractor EH.- There 
eh, n x h 2 0 was calculated numerically for a couple of SUSY models which resurrect higgsino 
m or wino m LSP and can yield acceptable 

Contrary to ref. ED, we focus on the range of the exponential-potential parameters, which 
ensures a late-time attractor together with an early KD regime (see sec. EH and can lead to a 
simultaneous satisfaction of several observational data (see secs. 12.41 and 14.211 . We then, present 
a “unified” (using the same independent variable) description of the cosmological evolution of 
the quintessence field and the x decoupling (see secs. 12.21 and 14.41) . The relevant equations are 
solved both numerically (see secs. 12.11 and 14.21) and semi-analytically (see secs. 12.21 and 14.41) 
and the results are compared with each other (see secs. 14. ll and 14 . 4l) . Finally, we demonstrate 
the crucial correlation between the D x /iq enhancement (w.r.t the one in the SC) and the 
quintessential density parameter at the eve of nucleosynthesis (see sec. IQ1) and we restrict 
the parameters imposing all the DE and CDM constraints (see sec. m without, however, to 
adopt a specific particle model. We showed that values of the quintessential density parameter 
at the former point close to its upper bound require the thermal-averaged cross section times 
the velocity of x to be almost three orders of magnitude larger than this needed in the SC. 

The framework of the quintessential cosmology is described in sec. El while our numerical 
and semi-analytical I\/iq calculations are displayed in sec. El Some numerical applications 
are presented in sec. EJ Finally, sec. El summarizes our results and discusses some open 
questions. Throughout the text and the formulas, brackets are used by applying disjunctive 
correspondence, natural units (h = c = Jzb = 1) are assumed, the subscript or superscript 0 is 
referred to present-day values and In [log] stands for logarithm with basis e [10]. 
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2. Quintessential Cosmology 

We briefly describe the equations which govern the evolution of the universe in the presence 
of quintessence (sec. ED, the phases which the quintessence field undergoes during its evolution 
(sec. ED and the requirements which a successful quintessential scenario is to satisfy (sec. ED - 

2.1. Relevant Equations 

According to the quintessential scenario, we assume the existence of a spatially 
homogeneous, scalar field q (not to be confused with the deceleration parameter 0) which 
obeys the Klein-Gordon equation. We below present its archetypal form and then we derive 
simplified forms which facilitate its numerical integration. Finally, we specify the used initial 
conditions and we define the useful extracted quantities. 


2.1.1. Initial form. The homogeneous Klein-Gordon equation in a cosmological set-up is 

(a) q + 3Hq + V q = 0, where (b) V = V 0 e~ Xq/m p (2.1) 

is the adopted potential for the field q, , q [dot] stands for derivative w.r.t q [the cosmic time, 
t] and H is the Hubble expansion parameter, 

(a) H = y/p q + p R + p u /VSrrip with (b) p q = i q 2 + V. (2.2) 

the energy density of q and m p = AIp/8n, where Mp = 1.22 x 10 19 GeV is the Planck mass. 
The energy density of radiation, p R , can be evaluated as a function of temperature, T, whilst 
the energy density of matter, p M , with reference to its present-day value: 

( a ) Pn = ^0 Sp * ^ and ^ PmR3 = P ° mR ° ( 2 ' 3 ) 
with R, the scale factor of the universe. Assuming no entropy production caused by the 
domination of another field (entropy production due to the q domination is not expected, since 
it does not couple to matter), the entropy density, s, satisfies the following two equations: 

2yr 2 

(a) sR 3 = Sp-Rp where (b) s = —g s * T 3 , (2.4) 

where subscript “p” represents a specific reference point at which the constant quantity sR 3 
is evaluated and g p *{T) [< 7 S *(T)] is the energy [entropy] effective number of degrees of freedom 
at temperature T. Their numerical values are evaluated by using the tables included in 
micrOMEGAs 1451 . originated from the DarkSUSY package |4]B] (recent improvements |47| do 
not affect essentially the results). 


2.1.2. Reformulation. The numerical integration of eq. ED is facilitated by converting the 
time derivatives to derivatives w.r.t the logarithmic time EHE3: 


t = In (R/Rq) = — ln(l + z) (=>- i: = H) 


(2.5) 


with 2 the redshift. Changing the differentiation, eq. ED turns out to be equivalent to the 
system of two first-order equations (prime denotes derivative w.r.t t): 


(a) Q = Hq' and (b) HQ 1 + 3 HQ + V q = 0 with (c) p q = -Q 2 + V. (2.6) 


In terms of t in eq. ED’ s and T can be expressed through the relations : 

.0 \ i/3 


(a) s = soe 3r and (b) T = Tq 


9s*\ 
9s* ) 


(2.7) 
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where eqs. E3t>) and E3t>) have been used. Similarly, p R and p M are elegantly cast in the 
form: 

/ o \ 4 /3 

( a ) PK = Pirfr(^r) e ” 4T and ( b ) PM = P° M e~ 3T - ( 2 - 8 ) 

i/p* \ys* j 

Eq. (12.8b ) was extracted by inserting eq. (ETTb ) in eq. (12.3b ). Eq. (12.8b ) can be derived by 
combining eq. (12.51) with eq. (IPb). 


2.1.3. Normalized form. An even more numerically “robust” .23. form of eq. (ESI can be 
achieved, if we introduce the following dimensionless quantities: 

( a ) ^m[r] = Pmib.]/Pc’ ( b ) Vo = V o /p 0 c and (c) q = q/V3m p . (2.9) 

Employing these quantities, eq. EU can be re-written as: 

(a) Q = Hq' and (b) HQ'+3HQ + V q = 0 with (c) H 2 = P q +P R +P M , (2.10) 
where the following quantities have been defined: 

(a) V = V 0 e~^, H = H/H 0l Q = Q/^~p ° and (b) p q = Q 2 /2 + E.(2.11) 
In our numerical calculation, we use the following values: 

p° = 8.099 x 10 ~ A7 hl GeV 4 and H 0 = 2.13 x 10 ~ A2 h 0 GeV, (2.12) 

with ho = 0.72. Also, = 0.29 and To = 2.35 x 10~ 13 GeV. Substituting the latter in 
eq. (12.3b ). we obtain p R = 8.04 x 10~ 5 . 

2.1.4. Extracted quantities. The solution of eqs. EZEHt O and Ennt>) allows us to calculate 
some measurable quantities which are used in order to test the quitessential model against 
observations (see sec. ESI). These are the density parameters of the g-field, radiation and 
matter 

G* = Pi/{Pq + p R + p M ) = Pi/H , where i = q, R and M, (2.13) 

respectively and the equation-of-state parameter (or barotropic index) of the (/-field, w q , 

w q = (q 2 / 2 - V)/(q 2 / 2 + V) = 1 - V/2p q . (2.14) 


2.1.5. Initial Conditions. In order to solve eq. E~mi) two initial conditions are to be specified: 
These could be the values of q and q 1 at an initial t, t : E3j. We take q( t : ) = 0 throughout our 
investigation, without any lose of generality. This is because possible use of q( t : ) = q l / 0 is 
equivalent as if we had q{ t : ) = 0 and rescaled V Q to V^exp (—\/3 XqJ 23) • This displacement 
influences just the choice of V Q determined from eq. ESS). 

On the other hand, the value of q'( t : ) is not a suitable initial condition for our purposes. 
This is, because we wish to focus on the regime (see eqs. (ITm (fTTTb ) and (OTTTTb )): 

flg(Ti) ~ ^ iQ ~ 1 ^ q'( Tj) ~ V2, (2.15) 

H i 

where we take p qi = p q (t, ) ~ Q 2 ( t : )/ 2 and H l = H( t : ). This means that Qi = Q(t : ) tends to 
infinity, since inserting eqs. (HTTTb ) and EB) into eq. Ennt) we can obtain: 


Q 



v + p k + p^ 

1 - q ' 2 /2 


( 2 . 16 ) 


In order to handle properly this subtlety, we find it convenient to define as initial condition, 
the square root of the kinetic-energy density of q at t I5 


= Qfa)/V2 ~ y/,p,~ ~ H r 


(2.17) 
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2.2. Quintessential Dynamics 

We can obtain a comprehensive and rather accurate approach of the q dynamics, following 
the arguments of ref. m- Namely, q undergoes the following three phases: 


2.2.1. Kination Dominated Phase. During this phase, the evolution of both the universe and 
q is dominated by the kinetic-energy density of q. Consequently, eq. EM) reads: 

(a) HQ' + 3HQ = 0 with (b) H = ^f q ~Q/V 2. (2.18) 

The former equation can be integrated trivially, with result: 

Q = 0ie -3(T_Tl) =>■ p q = p qi e~ 6 ^~ T ^. (2.19) 

Combining eq. (HTTKb l with nm we obtain: 

q~V 2(t—t : ) (=)> q = V 2) for t < t kr (2.20) 


with t kr , the point where the totally KD phase is terminated. This occurs, when: 


Pr/kr) Pij( t kr) 


~ + In 



( 2 . 21 ) 


where the right hand side of eq. (12^31) has been equated to the expression below at t = t kr 
( since t is close to t : we suppose that [g s *] does not vary from its value at z l . [//): 

/ 1 \ 4/3 

Pr = Paifr ( “ ) e' 4(T -b), with p RI = Pr (t t ). (2.22) 

£/p* Vi/s*/ 


Eq. (I2.4h l with reference point T t and eq. iU ) were employed in order to extract eq. (I2~22ll . 


2.2.2. Frozen-Field Dominated (FD) Phase. For t > t kr , the universe becomes RD but the 
evolution of q continues to be dominated by its kinetic energy density. Therefore, 

(a) HQ' + 3 HQ = 0 with (b) H = (2.23) 

Inserting eq. E2It>) into eq. GUBO and integrating the resulting one, we obtain: 

q = Qkr + P7u / 2 Q^ (Zkr - Zi) (1 - e" (T ' TK R)) for t kr < r < t fa (2.24) 

where g KR = q(t kr ) = In (p qi /p R1 ) /\/2 from eqs. (12.201) and (12.211) and t fa is specified in 
eq. (12^31) . It is obvious from eq. (I TM that q freezes at about t kf ~ t kr + 6 to the value: 

% - q KR + ^ (=> q = °) for ^kf < T < t fa- ( 2 - 25 ) 

Note that p q reaches its constant value, p qF = V(q F ), at t pl S> t kf such, that: 

<3 2 (t pl )/2 = V(q F ) => tt pl = Tj + Xq F /2V3-ln(y o /p qi )/6. (2.26) 

2.2.3. Attractor Dominated (AD) Phase. For t > t pl , p q becomes V dominated as in the case 
of inflation. Consequently, the evolution of q is described by the following: 

(a) 3HQ + V q ~0 with (b) H = ^ ~ (2.27) 

where p B is the dominant background-energy density of the universe with w B = 1/3 [0] for the 
RD [matter-dominated (MD)] era. As can be shown [32[, and has been extensively discussed 

E3E3G3II3DI, the system in eq. HU> admits: 
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(i) A global attractor for A > ^3(1 + w B ) with a fixed-point equation-of-state parameter 
w ip = w B and density parameter flq P = 3(1 + w B )/X 2 . This is the so-called self-tuning 
1231221132] case where the q evolution is insensitive to the choice of V Q . However, this 
case can be discarded (23 EU since, it fails to meet the observational data (see sec. EH). 

(ii) A late time attractor for A < y^3(l + w B ) with w q = A 2 /3 — 1 and = 1. As is pointed 
out 3 ESI an d we verify in sec. 14.21 the model can satisfy a number of observational 
constraints, for a reasonable set of initial conditions. 

Inserting eq. (HUS O into eq. tsza o and integrating the resulting equation, we obtain for the 
latter case: 

q = \z/V3 + \n(V o /p 0 q )/V?>X (=> q'= X/y/3) for t > t pa (2.28) 

where the transition from the FD to the AD phase occurs at the point t fa , which can be 
estimated by: 

t fa = v/6/A + V3q K JX - ln(V o /p 0 q )/X 2 . (2.29) 

The latter can be easily extracted by equating the values of the expressions in eqs. $TR 1) and 
(Egg)) for t = t fa . Employing eq. (15^51) . we can derive Q via eq. dmnt). Inserting it in 
the relation p q = Q 2 /(1 + w q ), which can be derived from eq. (12.141) . we arrive at the energy 
density of the late-time attractor: 

p A ~ with wf = A 2 /3 - 1. (2.30) 

2.3. Quintessential Requirements 

We briefly describe the various criteria that we impose on our quintessential model. 

2.3.1. KD “Constraint”. For the purposes of the present paper, we desire to focus our attention 
on the range of parameters which ensure an absolute [at least relative domination] of the q- 
kinetic energy at z 1 . This can be achieved, when: 

(a) Dq( tQ = 1 [(b) 0.5 < Q q ( tJ and (c) fl q (r t ) < 1] . (2-31) 

Ranges of parameters, which meet all the residual constraints of this section, not restricted by 
eq. are explored in ref. ESI- 

2.3.2. Nucleosynthesis (NS) Constraint. The presence of p q has not to spoil the successful 
predictions of Big Bang NS which commences at about t ns = —22.5 corresponding (see 
eq. (I2.7b ll to Tns = 1 MeV m ■ Taking into account the most up-to-date analysis of ref. m, 
we adopt a rather conservative upper bound on fl q (r NS ), less restrictive than that of ref. [49h 
Namely, we require: 

= n q (z NS ) < 0.21 (95% c.l.) (2.32) 

which corresponds to additional effective neutrinos species 5N U < 1.6 |48| . Let us note that 
extra contribution in the left hand side of eq. (1051) due to energy density of gravitational 
waves (GWs) created during the transition from the KD to RD era is negligible as we infer 
by explicitly applying the formulae of ref. ED- On the other hand, we do not consider 
contributions (potentially large J23) due to GWs generated during a possible former transition 
from inflation to KD epoch. The reason is that inflation could be driven by another field 
different to q and so, any additional constraint arisen from this period would be highly model 
dependent. 
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2.3.3. Coincidence Constraint. The present value of p q , ff ] q . must be compatible with the 
preferred range of eq. (int>). This can be achieved by adjusting the value of V 0 . Since, this 
value does not affect crucially our results (especially on the CDM abundance), we decide to 
fix p q to its central experimental value, demanding: 

Q° q = p° q = 0.73. (2.33) 

2.3.4. Acceleration Constraint. A successful quintessential scenario has to account for the 
present-day acceleration of the universe, i.e. [1_, 

- 1 < w q {0) < -0.78 (95% c.l.). (2.34) 

In addition, since the string theory disfavors the eternal acceleration, it would be desirable to 
demand > —1/3 TOT . However, in the case of the used potential, we did not succeed to 
achieve compatibility of the latter optional restriction with eq. dunn) , in accordance with the 
findings of ref. J2fj). 

2.3.5. Residual Constraints. In our scanning, finally, we take into account the following less 
restrictive but also non-rigorous bounds, which, however, do not affect crucially our results: 

(a) -50 <t : <- 36 and (b) TTi < 10 56 . (2.35) 

The lower bound of eq. (12.35b ) comes from the gravitino constraint EH which provides an 
upper bound on the reheat temperature, Trh < (10 9 — 10 10 ) GeV. This can be translated to a 
lower bound on t : , through eq. EZb). However, this bound may not be so reliable, since there 
is no thorough investigation of the gravitino constraint within the context of quintessential 
cosmology, to date. Also, since we do not study the evolution of the universe before the 
commencement of the KD era, we wish to liberate our O x /ig calculation from this ignorance. 
To this end, we demand t : to be lower than the upper bound of eq. (I2.35b l. This corresponds 
to the onset of the Boltzmann suppression of the y-nurnber density (see sec. EH for mass of y 
equal to 500 GeV (see eq. (EH)). The bound of eq. (I2.35b l comes from the COBE constraints 
M. on the spectrum of GWs produced at the end of inflation m- 

3. CDM Abundance in the Presence of the KD Phase 

We assume that the CDM candidate, y, maintains kinetic and chemical equilibrium 
(see below) with plasma, is produced through thermal scatterings and decouples (being non- 
relativistic) during the KD epoch. Our theoretical analysis is presented in sec. 13.11 and its 
numerical treatment in sec. 13.21 Useful approximated expressions are derived in sec. 13.31 

3.1. The Boltzmann Equation 

Since the y particles are in kinetic equilibrium with the cosmic fluid, their number density, 
n x , satisfies the following Boltzmann equation: 

h x + 3 Hn x + {av) (ra 2 — n® q2 ) = 0, (3.1) 

where H is given by eq. EHO, {av} is the thermal-averaged cross section of y particles 
times the velocity and n x l is the equilibrium number density of y, which obeys the Maxwell- 
Boltzmann statistics: 

nl q (x ) = -—t 3 / 2 e~ 1 / x P 2 (l/x), where x = T/m x 
A (27r)^/^ a 


(3.2) 
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with m x the mass of x- We pose g = 2 for the number of degrees of freedom of x and 
P n (z ) = 1 + (4n 2 —1)/8 z is obtained by asymptotically expanding the modified Bessel function 
of the second kind of order n. Note that non-chemical-equilibrium production of x’s requires 
(crv) < 1(T 20 GeV -2 . Since such a value is well below the usually obtainable values 
HU 00 EH], we do not consider further this possibility, here. 


3.2. Numerical Calculation 

Following the strategy of sec. 12.1.31 we introduce the dimensionless quantities: 

pleq] = p [eq]/ p 0' w here pN = m x n ^. (3.3) 

In terms of these, eq. m takes the following master, for numerical manipulations, form : 

Hp x + 3Hp x + (crv) (p x — pf) /rh x = 0, (3.4) 

where H is given by eq. (Hunt) and the following quantities have been defined: 

(crv) = \ffJl(av) and fh x = m x /V3m p . (3-5) 

Eq. (13.41) can be solved numerically with initial condition p x (z B ) = / 0 ^ q (r B ), where t b 
corresponds (see eq. COM) to the beginning (x = 1) of the Boltzmann suppression of p f 
Since t : < t b , the integration of eq. 03) is realized from t b down to 0. Finally, Q x Kq can be 
easily found, via the relation: 

= P x m l (3.6) 


3.3. Semi-Analytical Calculation 

The aim of this section is the calculation of O x /iq based on the already obtained semi- 
analytical expressions of sec. 12.21 The procedure is described step-by-step below. 


3.3.1. Reformulation of the Boltzmann Equation. Introducing the variables yt eq l = njf/ s 
0 ESI (in order to absorb the dilution term) and converting the derivatives w.r.t t, to 
derivatives w.r.t t, eq. m can be rewritten as: 

HY' = -(av) ( Y 2 - T eq2 ) s, (3.7) 

where eq. E30 has been also utilized. Substituting eqs. and 0H in eq. 02t) and 

ignoring the negligible, during the x decoupling, contribution of p M . H can be expressed as: 


Equivalently for t, 
we obtain: 


(a) H = \J p R g q /3mp, where (b) g q = 1 + r q with (3.8) 

1 / n 4/3 

r q = r,— ( ) e _2 ( T_T i) an d r, = (3.9) 

9 P * \9sJ Pm 

> taking as reference point t ns instead in eqs. 021 and (I2T91) . 

fcV /3 


NS 

tfp* 


9p* \g^ s J 


; 2 ( T t ns) with 


r NS NS 

' R 


of 


1-Of 


(3.10) 


and the superscript NS denotes the values of the several quantities at t ns . Inserting eqs. 01 , 
030 and 03) into eq. 03 this can be cast in the following final form: 

Y' = yJ?™L (; Y 2 - y eq2 ) , where (3.11) 

\[9q 

( 0 \ ^ . 2 

t) (f) with Hl= 


(3.12) 
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3.3.2. The freeze-out procedure. In the case of the equilibrium x production, an accurately 
approximate solution of eq. cmi) can be achieved, introducing the notion of freeze-out 
temperature, T F = T( z F ) = x F m x [9, [12], which allows us to study eq. (13.111) in the two 
extreme regimes: 

• At very early times, when z <C t f , \’s are very close to equilibrium. So, it is more 
convenient to rewrite eq. (ETTTl) in terms of the variable A(t) = Y(z) — Y eq (z) as follows: 

A' = -Y eq/ + y (av) A (A + 2Y eq ) /^g~ q . (3.13) 

The freeze-out point t f can be defined by 

A(t f ) = d F y eq (T F ) =► A(T F )(A(T F )+2y e< i(T F )) = Sf(5 f +2) Y eq2 (z F ), (3.14) 

where d F is a constant of order one, determined by comparing the exact numerical solution of 
eq. (THTl) with the approximate under consideration one. Inserting eqs. m into eq. dim 
we obtain the following equation, which can be solved w.r.t t f iteratively: 

(lnY eq ) (t f ) = y F (av)5 F (5 F + 2)Y eq (z F ) / ^(6 F + 1) with (3.15) 

y,=y(t,) and (l„r«)V)=*'(T-|--|i + i AL- y ). (3.16) 

• At late times, when t> t f , h > Y eq and so, Y 2 — Y eq2 ~ Y 2 . Inserting this into 
eq. cmi) and integrating the resulting equation from t f down to 0, we arrive at: 

(a) Yq = (YpT 1 + J F ) 1 , where (b) ,/ F = f dr — ^° V ^ and (3-17) 

Jo y/dq 

Yp = (Jp + 1) y eq ( Tp ) with Y eq (z) = — e ~ l/x -^(I/t), (3.18) 

9s* y 2 

where the x — t dependence can be derived from eq. (1277b T Although not crucial, a choice 
d F = 1.2 =F 0.2 assists us to approach better the precise numerical solution of eq. nm . 

3.3.3. The CDM abundance. Our final aim is the calculation of the current x relic density, 
which is based on the well known formula m- 

Q x = p° x /p° c = m x s 0 Y 0 /p° c => Q x h 2 0 = 2.741 x 10 8 Y 0 m x / GeV. (3.19) 

The presence of g q > 1 in eq. (13.151) and, mainly, in eq. (ITTTTTb ) reduces J F , thereby increasing 
the value O x /ig w.r.t the one obtained in the SC (i.e. with g q = 1), O x /ig| gc . The resulting 
enhancement can be estimated, by defining the quantity HQ: 

Afl x = (O x / lo -^^olso)/^^ol s c- (3.20) 

3.3.4. The variation of AO v . The variation of AO x w.r.t the free parameters can be designed 
by simplifying the formulas above. In particular, z F and Ad x can be roughly estimated as: 

(a) z F ~ —2ln( x /m x /(av)) and (b) AO x ~ J F | SC / Jf ~ 1 ~ yfdq ~ 1 ~ (3.21) 

where we kept only the most important terms in eqs. cm, mB and (HDH) . Also, we 
have taken into account eqs. (I3.17b h from which we extracted J F ~ (av)e~ r F /T F | SC ~ 
(ov)e~ T v (for constant (av) = a). 

Armed with these formulas, we can explain that Ally increases as: (i) g q increases (for 
fixed m x and (av))-, this is obvious from eq. (I3.21b l. (ii) (av) decreases (for fixed m x and g q ); 
indeed, from eq. (EHTb T decrease of (av) results to a decrease of t f , which in turn, causes an 
increase of Afl x . (iii) m x increases (for fixed (av) and g q ); indeed, as shown from eq. (I3.21b l. 
an increase of m x generates a decrease of t f which increases Afl x . 
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4. Applications 

Our numerical investigation depends on the parameters: 

A, t i; H lt m x , (av). 

For ease of reference we call the three first parameters q parameters, whereas the two later, 
CDM parameters. Recall that we use q{ tJ = 0 throughout and V a is adjusted so that eq. (12.4411 
is satisfied. Nonetheless, for definiteness and clarity, we give the used value of V Q in the explicit 
examples of figs. 3 Eland01 hr general, V 0 ranges between about 1 and 10 33 , increases with 
A or H l and turns out to be (t : ,^^- independent, for fixed A and fi q (Tj) < 1. 

As regards the CDM parameters, we have to clarify that (av) can be derived from m x 
and the residual (s)-particle spectrum, once a specific theory has been adopted. However, to 
keep our presentation as general as possible, we decide to treat m x and (av) as unrelated input 
parameters (following our strategy in ref. [HjaJ). Specifically, keeping in mind that the most 
promising CDM particle is the LSP, we focus our attention on the range: 

200 GeV < m x < 500 GeV. (4.1) 

Taking into account the experimental constraints on the SUSY spectra of several SUSY models 
(see, e.g., fig. 23 of ref. [3]), we adopt a rather restrictive lower bound on m x which, however, 
ensures us that the range of eq. 63 is valid even in the most constrained cases. The upper 
bound in eq. 63 is imposed in order the analyzed range to be possibly detectable in the 
future experiments (see, e.g. ref. |H3i)- On the other hand, we isolate the two extreme cases 
which we encounter when we use the non-relativistic expansion in order to calculate (av) (the 
method gives, in general, accurate results far enough from s-poles and thresholds mmy- 

(a) (av) = a or (b) (av) = bx. (4.2) 

The x dependence in eq. 63)) emerges in the case of a bino LSP m without coannihilations 
(CANs), whereas eq. (14.2b 1 is extracted in the majority of the residual cases 0 - 0 , m-m- 
The values of a and b can be restricted by applying the bounds of eq. 63- Comments on the 
naturalness of the required values are given in sec. 14.4.11 

The presentation of our results begins with the description of the evolution of the various 
quintessential quantities in sec. 14.11 In sec. 14.21 we present the ranges of the quintessential 
parameters, allowed by the constraints of sec. 12.41 In sec. roi we investigate the behaviour of 
the f\h.Q enhancement and finally, in sec. 14.41 we present areas compatible with eqs. 63- 

4.1. Evolution of the Quintessential Quantities 

We illustrate the evolution of the various quintessential quantities presenting diagrams 
where in the x-axis, t= — ln(l + z) varies from T t down to late times, e.g. 10 122 G23> EU 

In fig. 3( a ), we display q and q' versus t for Xj = —41, H l = 10 40 and A = 0.5 ( V Q = 
1.6 x 10 8 ). Solid lines [crosses] are obtained by numerically solving eq. (12.101) [applying the 
analytical expressions of sec. 63 Despite their simplicity, our semi-analytical expressions in 
eqs. OB . 63i and 633), reproduce impressively the numerical evolution of q and q'. 

In fig. 3(b), we draw log ^ versus t for t : = —41, H l = 10 40 and two “extreme” (see 
figs.|31-(a) and (b)) values of A, A = 0.1 ( V a = 33.5) (solid lines) or A = 1.15 (V a = 1.35 x 10 19 ) 
(dashed lines). For i = q (bold black lines), we show log p q , computed by inserting in eq. (I2.11b l 
the numerical solution of eqs. 6330 and (I2.10b l. For i = A (thin black lines), we show log/^ 
derived from eq. 633) - For i = R + M (light grey line), we show log p R+M , which is the 
logarithm of the sum of the contributions given by eqs. 630 and 63). 
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Figure 1: The evolution as a function of r for % l = —41 and H I = 10 40 of the quantities: 
q and q' for X = 0.5 ( V 0 = 1.6 x 10 8 ) (crosses are obtained by our analytical expressions) 
(a), log with i = q (bold black lines), A (thin black lines) and R+M (light grey line) for 
X = 0.1 (Vo = 33.5) (solid lines) or X = 1.15 (Vo = 1-35 x 10 19 ) (dashed lines) (b), w q 
(dark gray line) and fli with i = q (black line), R (light gray line) and M (gray line) for 
X = 0.5 (Vo = 1.6 x 10 s ) (c). 

In fig. row, we plot w q (dark gray line) and fij with i = q (black line), R (light gray line) 
and M (gray line) versus t, for t [ = —41, H 1 = 10 40 and A = 0.5 (V 0 = 1.6 x 10 8 ). The y -axis 
quantities are computed by inserting in eqs. dznp and (TTHTl) the numerical values obtained 
by eqs. C3It), (THTb l and CEB) . 

Analyzing comparatively figs. EKa), (b) and (c), we can demonstrate the characteristic 
features of the cosmological history in the presence of q. In particular: 

(i) For T : < t < t kr = —26.2, the universe undergoes the KD era. The field q increases 
according to eq. TFM along the left inclined part of the curve in fig. EKa) and p q decreases, 
more steeply than p R , according to eq. (Em along the left inclined part of the black solid 
curve in fig. nKb). During this period, w q = 1, as shown in fig. ^-(c). This era terminates 
at t kr , where an intersection of p q [f2 g ] with p R [Dr] is observed in fig. H-(b) [fig. EK C )]- 
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Figure 2: The evolution of the quantities log p i with i = q (black lines ) and R+M (light 
grey line ) as a function of r for X = 0.5 and : Hi = 10 40 and r 1 = —40.025 ( V 0 = 1.6 x 10 9 ) 
/t, = —48.64 (Vo = 2 . 3 )] (solid [ dashed ] line ) (a), t, = —40 and H 1 = 9.3xl0 39 ( V 0 = 1.6xl0 9 ) 
[Hi = 3.2 x 10 32 (V 0 = 2 . 3 )] (solid [ dashed ] line ) (b). 

(ii) For t kr = —26.2 < t < t fa = 0.67, the universe undergoes successively the RD era and 
then the MD era until the re-appearance of DE. More precisely, q freezes to its constant 
value q F — 22.4 according to eq. (IQTli (or eq. (IQKll for t > t kf = —20.2) along the 
horizontal part of the curve in fig. ID(a). V fig increases towards 1 along the light 
gray line of fig. □(c) and then decreases until fl M = fl R at r eq = —8.16, where a slight 
kink is observed on the light gray line of fig. HK b). F or t kf <t< t pl = -10.22, p q < p R 
continues to decrease steeply according to eq. (J2H, along the left inclined part of the 
black solid curve in fig. □(b), while w q continues to be 1 as shown in fig. □ (c). On the 
other hand, for t pl < t < t fa , p q freezes to its constant value p qF = V(q F ) = 0.62 in 
fig. □(b), while w q transits from 1 to -1 as shown in fig. ^( 0 ). At present, we obtain 
w q ( 0) = —0.96, within the limits of eq. (12.441) . 

(iii) For t > t fa = 0.67, the universe undergoes a (/-dominated phase. The field q increases 

according to eq. along the right inclined part of the curve in fig. □(a), p q decreases 

according to eq. EDI along the right inclined parts of the curves in fig. □ (b) for A = 0.1 
and 1.15, while w q in fig. □ (c) tends to its fixed-point value for A = 0.5, = —0.92. As 

shown in the same figure flg P = 1, which identifies the late-time attractor. 

Note that, contrary to the case with A > 2 @U, the variation of A does not affect essentially 
the position of the FD plateau but only changes the inclination of the AD curve. 

The dependence of the p q evolution on t y [H i ]. can be easily concluded from fig. □(a) 
[fig. 13(b)], where we plot for A = 0.5, log/jg and log/0 R+M versus t. In fig. 13(a) [fig. □(b)] we 
take Hi = 10 4 ° and t : = —40.025 (V Q = 1.6 x 10 9 ) (solid line) or = —48.64 (V 0 = 2.3) (dashed 
line) [t : = —40 and H 1 = 9.3 x 10 39 (V 0 = 1.6 x 10 9 ) (solid line) or H\ = 3.2 x 10 32 (V Q = 2.3) 
(dashed line)]. It is obvious that increasing t : or H l , the left, black inclined line of the q-KD 
regime moves to the right and consequently, both Q q ( t,) and Dg(T NS ) increase. So, an upper 
[lower] bound on t : or H l can be extracted from eq. (IOH [eq. (1011) 1 (see fig. ®. The 
saturation of these inequalities is the origin of the chosen lower [upper] or H l in fig. [3 
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Figure 3: The allowed (shaded) areas on the t, — A plane for H 1 = 10 40 (a), log H 1 — A plane 
for Tj = —40 or t 1 = —50 (b) and t 1 — logH 1 plane for A = 0.5 (c). Lined are, also, the areas 
allowed by eq. 1 1. 21) for m x = 350 GeV and (ov) = 10 -7 GeV -2 . 

4.2. Imposing the Quintessential Requirements 

We proceed, now in the delineation of the parameter space of our quintessential model. 
Agreement with eq. dUSH) entails 0 < A < 1.15 (see also ref. where less restrictive upper 
bound on w q ( 0) was imposed). This range is independent on t : and H l as is shown in figs. 01(a) 
and 0Kb). In these, we depict respectively the allowed (shaded) regions on the t : — A plane 
for H l = 10 40 and on the logH 1 — A plane for = —40 or = —50. In fig. OKc), we design 
the allowed area on the t : — logl/j plane. Although this plot is constructed for A = 0.5, it 
is obviously A independent. The dark [light] shaded areas fulfill eq. (I2.31b l [eq. (12.31b ) and 
(12.31b )]. The right [left] boundaries of the allowed regions in figs. OKa) and 0Kb) are derived 
from eq. (12.321) [eq. (12.31b )]. The same origin has the upper [lower] boundary of the allowed 
region in fig. EKO, whereas the left and right boundaries come from eq. (12.35b ). So, for a 
reasonable set of (Athe exponential quintessential model can become consistent with 
the observational data nacsj. The construction of the ruled areas is explained in sec. 14.4.21 
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(a) 



x = - ln(1+z) 

(b) 


Input Parameters 

A 



V 0 

m x (GeV) 

(<ru)(GeV 2 ) 

0.5 

-41 

10 4 ° 

1.6 x 10 8 

350 

2.92 x 10“ 7 

Output Parameters 

^q( T ns) 


t f Isc 

t kr 

n x h 2 0 

AD X 

0.0007 

-32.9 

-32.6 

-25.7 

0.115 

126.8 


Figure 4: The evolution as a function of t of the quantities p^/s (dotted line) and p x /s 
(bold [thin] solid lines) for the quintessential scenario [SC] (a) and log p i with i = q [i = R\ 
(black [light grey] line) (b) for the input quantities listed in the table above. 


4.3. The D x /ig Enhancement 

The investigation of the D x /ig enhancement is the aim of this section. In sec. 14.3.11 we 
illustrate the x decoupling during the KD epoch and in sec. 14.3.51 we examine the dependence 
of the f\/tQ increase on D,j(t ns ). Finally, in sec. 14.3.31 we compare the results of our numerical 
and semi-analytical Q x /iq calculations. 

4.3.1. The x decoupling. The x decoupling during the KD era is instructively displayed in 
%s.[J(a) and (b). In fig. [1(a) we depict /0 x q V s (dotted lines) and p x /s [p x /s| sc (for g q = 1)] 
(bold [thin] solid lines) versus t. In fig. BKb), we plot logp q [logp R ] (black [light gray] line) 
versus t. The needed for our calculation inputs and some key-outputs are listed in the relevant 
table. For better comparison, we give, also, the point of the % decoupling, in the case of the 
SC (g q = 1), t f | sc . In the present case, the x decoupling is realized deeply within the KD 
regime, t : < t f < t kr and t f < t f | sc . By adjusting ( av) we extract the central D x hg in 
eq. O- The presence of the KD era causes an efficient D x /ig enhancement, AD X = 126.8 
(fl x /iQ| sc = 0.0009). Note that the condition t f < t kr is indispensable in order to obtain 
sizable A Q x . This can be understood as follows. From eqs. Tn n and (THUl) we obtain 
t kr ~ t ns + ln r Ng /2. If we demand t kr < t f , we obtain D g (T NS ) < 7.7 x 10 -1 °, which causes 
a very weak AD X < 2.2. Finally, the phenomenon of re-annihilation m is not observed in 
this context. This is, because in our case H smoothly evolves from its KD to RD behaviour - 
see eq. (ESI) - and does not sharply drop after the x-decoupling as i n the case of ref. E3. 
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Af2 x 


m x = 350 GeV, 


(av) (GeV -2 ) 


10 -10 

10 -7 

0.21 

2242 

1769 

0.001 

187.7 

150.4 

10 -7 

3.42 

2.62 

10 -11 

0.0036 

0.02 


Figure 5: Lines with constant S7 9 (t ns ), indicated in the plot, in the allowed t 1 — \ogH 1 plane. 
The corresponding constant values of Afor m x = 350 GeV and (ov) = 10 -10 GeV -2 or 
(crv) = 10 -7 GeV -2 are also listed in the table. 


4.3.2. The dependence of Afl x on Q. q (z NS ). As is shown in figs. 121(a) and[2-(b), the position 
of the inclined left part of the black line (corresponding to log p q ) is affected crucially by a 
possible variation of z l or H 1 but not of A. Therefore, £l q (z^ s ) and consequently, A12 x (see 
eqs. EH> and (13.10H 1 depend on z x or H 1 but not on A (contrary to the case of ref. .HIT 
Moreover, the dependence of A12 x on z x or H 1 can be expressed exclusively as a single-valued 
function of Q q (z NS ), since only g q is involved in the O x /iq calculation (see eqs. (13.191) and 
(13.171) 1. This is illustrated in the fig. El where we depict iso-fig( t ns ) lines on the allowed 
region of the z 1 — log H l plane, presented in fig. 01(c). Along these lines, A Q x remains, also, 
constant (indicated in the table of fig. EJ for fixed m x = 350 GeV and {crv) = 10 10 GeV 2 
or (crv) = 10“' GeV -2 . Consequently, the number of the free q parameters (t I5 H x ), which 
determine Aff x , can be reduced by one and replaced by Q q (z NS ). 

4.3.3. Numerical Versus Semi-Analytical Results. The validity of our semi-analytical approach 
can be tested by comparing its results for Aff x with those obtained by the numerical solution of 
eqs. (E?~l . In addition, useful conclusions can be inferred for the behavior of Aff x as a function 
of our free parameters, m x , (crv) and flg(% s ). Our results are presented in fig. El The solid and 
dashed lines are drawn from our numerical code, whereas crosses are obtained by employing the 
formulas of sec. 13.31 with <5 f = 1-28 [<5f = 1-35] for (crv) = a [(crv) = bx]. In figs.El( a i) EK a 2 )] 5 
we present AQ X versus ^(t^) for (crv) = 10 -10 GeV -2 [(crv) = 10 -10 x GeV -2 ]. We take 
m x = 200 GeV [m x = 500 GeV] (light [normal] grey lines and crosses). In fig. El(b), we plot 
AQ X versus a = (crv) [a = (crv)/x] for m x = 350 GeV and (crv) = a [(or) = ax] (solid [dashed] 
lines). In fig.EK c ) 5 we depict AQ X versus m x for (crv) = 10 -10 GeV -2 [(crv) = 10 -10 x GeV -2 ] 
(solid [dashed] lines). In the last two cases, we use £l q (z NS ) = 0.001 [n g (T NS ) = 0.1] (light 
[normal] grey lines and crosses). As we anticipated in sec. 13.3.41 Af2 x increases when Q q (z NS ) 
(see figs. El( a i) and (a 2 )) or m x (see fig. El( c )) increases and when (crv) decreases (see fig. El 
(b)). From figs. El(b) and El( c ) is, also, deduced that AI2 X increases more drastically in the 
(crv) = bx case than in the (crv) = a case for a = b and fixed m x and Q q (z^ s ). Evident is, 
finally, the agreement between numerical and semi-analytical results. 
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(b) (c) 

Figure 6: Af2 x versus 12 9 (t ns ) for m x = 200 GeV [m x = 500 GeV] (light [normal] grey lines 
and crosses) and (ov) = 10 -10 GeV -2 [(era) = lO -10 * GeV -2 ] (ai [a 2 ]). Also, Afi x versus 
a = (ov) [a = (av)/x] for m x = 350 GeV and (ov) = a [(err) = ax] (solid [dashed] lines) (b) 
and Af2 x versus m x for {ov) = 10 -10 GeV -2 \{ov) = 10 -10 * GeV -2 ] (solid [dashed] lines) 
(c). We take S2 9 (t ns ) = 0.001 [£2 9 (t ns ) = 0.1] (light [normal] grey lines and crosses). The 
solid or dashed lines [crosses] are obtained by our numerical code [semi-analytical expressions]. 


4.4. Imposing the CDM Requirement 

Requiring to be confined in the cosmologically allowed range of eq. m , one can 

restrict not only the CDM parameters (see subsec. 14.4.Ill but also the q parameters, A, z l 
and H l (see subsec. E33 or Dq(r NS ) (see subsec. EH;. The data is derived exclusively by 
the numerical program. Let us note, in passing, that bounds arisen from eq. Ob) , are more 
rigorous than those originated from eq. Qb), since other production mechanisms of y’s may 
be activated (SB, m and/or other CDM candidates jjj may contribute to Dcdm- 

4.4.1. Constraining the CDM parameters. Fixing the q parameters, we can derive restrictions 
on the CDM parameters. Namely, in fig. 0we construct the allowed regions on the m x — (av) 
plane for (av) = a (a, ai, 32 ) or on the m x — (av)/x plane for (av) = bx (b, bi, b 2 ). As we 
showed in subsec. 14.4.21 the q parameters can be replaced by D q (T NS ). So, in figs. 0-(ai) and 





















<gv> = a (x 10 6 GeV' 2 ) <gv> = a (x 10 7 GeV' 2 ) <gv> = a (x 10 9 GeV' 2 ) 
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Figure 7: The allowed region on the m x — (av) plane for (av) = a (a, ai, a 2 ) or on the 
m x — (av)/x plane for (av) = bx (b, bi, b 2 ). We take fi q (T NS ) = 0.001 [D 9 (t ns ) = 0.1] in fig. 
(ai, bi) [fig. (a 2 , b 2 )7 whereas, for the sake of comparison, we consider the SC (H 9 (t ns ) = 0) 
in figs, (a) and (b). 
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(a) 


(b) 


Figure 8: The allowed regions for m x = 200 GeV [m x = 500 GeV] (light [normal] grey areas) 
on the (av) — IUUns) plane for (av) = a (a) or on the (av)/x — Q. q (z NS ) plane for (av) = bx 

(b). 


(bi) [figs. 0- (a 2 ) and (b 2 )] , we fix Q q (r NS ) = 0.001 [H 9 (t ns ) = 0.1], whereas in figs. 0- (a) and 
(b), we consider, for better reference, the SC, with D g (T NS ) = 0. The upper [lower] boundaries 
of the allowed areas are derived from eq. (II. 2b ) [eq. (11.2b )]. This is due to the fact that f\/ig 
is inverse proportional to (av) as is obvious from eqs. mm and (1X171) and so, decreases 
as (av) increases (contrary to the case of the non-equilibrium x production (35, [3BJ). 

We observe that with H 9 (t ns ) = 0.001 [fl g (T NS ) = 0.1], agreement with eq. (11.21) entails 
almost two [three] orders of magnitude higher (cru)’s than those required in the SC. Also, due 
to the presence of Q q (T NS ) > 0, D x //q increases with m x more dramatically than in the case of 
the SC (fi q (T NS ) = 0), illustrated in figs. 0(a ) and (b). This effect is more straightened in the 
(av) = bx case, as is seen in figs. 0(bi) and (b 2 ). 

The requisite high values for (av) (almost unnatural in the (av) = bx case) can be obtained 
by resorting to SUSY models which ensure A-pole effects £DE2 or “gaugino-inspired” CANs 
157 , 5S1 I59| . as in the applications |43l I44j of ref. |41| . Less efficient augmentation of (av) 
can be achieved by lowering the masses of the CDM candidates in ED models 001 or 
by employing sfermionic CANs 155 .. Eli 00 in SUSY models. Consequently, the constrained 
minimal SUSY model m , although tightly restricted even in the SC mm, can become 
consistent with a quintessential KD period, e.g. applying the A-pole effects EHE2- 

4.4.2. Constraining further the q parameters. Fixing the CDM parameters to naturally 
obtainable values, m x = 350 GeV and (av) = 10 -7 GeV~ 2 (which yield Q x Iiq\ sc ~ 0.0025), 
we can constrain further the q parameters, which are already constrained by the quintessential 
requirements in subsec. 14.21 The regions consistent with the achievement of eq. o are ruled 
in fig. 01 As expected from the argument of subsec. 14.3.21 D,j(t ns ) turns out to be constant and 
equal to 8.5 x 10 -5 [3.6 x 10 -5 ] along the right [left] boundaries of the ruled areas in fig. El-(a), 
fig. 0(b) and along the inclined upper [lower] boundary of the ruled area in fig. 0(c). If we 
had imposed only the bound from eq. (Obh we would have obtained obviously much wider 
allowed regions bounded from the upper boundary of ruled area and the lower boundary of 
the light shaded area. 















19 


Quintessential Kination and CDM Abundance 


4.4.3. Constraining further D q (T NS ). Since H x /ig depends exclusively on D q (T NS )’s for fixed 
m x and ( av }, it would be interesting to delineate the allowed parameter space on the 
(av } — flq(T NS ) \(av)/x — D q (T NS )] plane for (av) = a [(av) = bx\ and fixed m x . This 
aim is realized in fig. [HI- (a) [fig. EKb)]. The light [normal] grey regions are constructed for 
m x = 200 GeV \m x = 500 GeV]. Lower m x s require lower (av)' s, since Af2 x decreases 
with m x as we explain in sec. 14., 3., 31 Also, the upper [lower] boundaries of the allowed 
areas are derived from eq. (II ,2h l [eq. (II .2h l]. for the reason already mentioned in sec. 14.4.1 1 
Consequently, when the CDM parameters are given, restrictions on Q q (z NS ) supplementary to 
those from eq. can be derived from eq. (E2J. 

5. Conclusions-Open Issues 

We studied the cosmological evolution of a scalar field q which rolls down its exponential 
potential ensuring an early KD epoch and acting as quintessence today. We then investigated 
the decoupling of a CDM candidate, %, during the KD epoch and calculated D x /ig. We 
solved the problem (i) numerically, integrating the differential equations which govern the 
cosmological evolution of q and the %-nurnber density (ii) semi-analytically, producing 
approximate relations for the former quantities. The second way facilitates the understanding 
of the problem and gives, in all cases, accurate results. 

The parameters of the quintessential model (A, t : , were confined so as 0.5 < flq(Tj) < 1 
and were constrained by using current observational data originating from nucleosynthesis, the 
acceleration of the universe and the DE density parameter. We found 0 < A < 1.15 and that 
there are reasonably allowed regions on the (t : , L/Q-parameter space. We also showed that 
D x /iq increases w.r.t its value in the SC with fixed m x and (av). We analyzed the variation of 
this enhancement, AD X , w.r.t (\,r 1 ,H 1 ) demonstrating that it can be expressed as a function 
of D q (T NS ). We, also, found that Af l x increases with Q q (t ns ) and m x and as (av) decreases. 
It is, also, larger in the (av) = bx case than in the (av) = a case for a = b and fixed m x and 
D q (T NS ). By enforcing the CDM constraint, D q (T NS ) close to its upper bound requires almost 
three orders of magnitude larger (av)' s than those required in the SC for fixed m x . 

Our formalism could become applicable to other more elaborated quintessential models 
HHUSHEZI Also, it could be easily extended, in order to include coannihilations and/or A-pole 
effects for the fl x /ig calculation in the context of specific SUSY or ED models. In the latter case, 
novel deviations m from the SC arise, which could be similarly analyzed (although the brane- 
tension is to be rather low in order numerically visible changes on the fQ/ig calculation to be 
observable EH)- Also, low-reheating scenaria •E MI IE] could become extremely appealing 
in the presence of quintessence, since they succeed to reduce I2 x /ig without need of tuning 
the particle-model parameters (their coexistence with the quintessential evolution deserves 
certainly deeper investigation ESI). On the other hand, the D x /ig enhancement is welcome for 
wino or higgsino LSPs 

jHEHlEEi, which yield D x /ig lower than the expectations in the SC, 
and so, Afl x can drive fl x /jg to the correct value m ■ In the same time, relatively high direct 
detection rates can be produced without invoking the questionable normalization ESI of the 
proton-nucleus cross section. 
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